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Dark energy and dark matter are the dominant sources in the evolution of the late universe. 
They are currently only indirectly detected via their gravitational effects, and there could be a 
coupling between them without violating observational constraints. We investigate the background 
dynamics when dark energy is modelled as exponential quintessence, and is coupled to dark matter 
via simple models of energy exchange. We introduce a new form of dark sector coupling, which leads 
to a more complicated dynamical phase space and has a better physical motivation than previous 
• mathematically similar couplings. 
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Ch ■ I. INTRODUCTION 



Observations are providing increasingly compelling evidence that the expansion of the Universe is accelerating, 
^sq' . driven by "dark energy" (see e.g. [l[ for recent results). The simplest model of dark energy is a cosmological constant 
A, representing the vacuum energy density, and this model provides a very good fit to a range of independent 
observations. However, there is no satisfactory theoretical explanation for the very small value of A. Furthermore, 
the A model suffers from a fine-tuning, or "coincidence" , problem - why is the dark matter density comparable to the 



I , vacuum energy density now, given that their time evolution is so different? 



5-H . 



If the dark energy evolves with time, this may alleviate the coincidence problem. The simplest models of evolving 
dark energy are light scalar fields, known as "quintessence" . If the quintessence is coupled to the dark matter, then 
this may be able to account for the similar energy densities in the dark sector today. A decisive way of achieving 
similar energy densities is if the coupling leads to an accelerated scaling attractor solution, with 

\0\ »dark energy ^ ^^^^ ^^^^ ^ ^ 

^ 'dark matter 



in 



o 

oo 
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In this case, the coincidence problem is reduced to a simple choice of parameters to match fldark energy /^^dark matter 
to observations. Since the accelerated scaling solution is an attractor, no fine-tuning of initial conditions is needed. 
However, the dynamics that produces such scaling in the dark sector may have other undesirable consequences. 
Here we study quintessence with an exponential potential, 

V{lp) = Vo exp {-nXip) , := SttG , (2) 



^ . , where A is dimensionless and Vq > 0. The dynamics of a universe with exponential quintessence and an uncoupled 
■ perfect fluid have been studied, and these models do not admit late-time accelerated scaling attractors that satisfy 
Eq. ([1]) When we introduce a coupling between the quintessence and dark matter, accelerated scaling attractors 
are possible However, in some models, this is achieved at the expense of introducing other problems which can 

rule out the model Q. 

A general coupling between a quintessence field ip and dark matter (with density pc) may be described in the 
background by the balance equations, 

Pc = -iHpc - Q , (3) 
p^ = -3H{l + w^)p^ + Q . (4) 
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Here Q is the rate of energy density exchange in the dark sector, and 

Q 

The dark energy equation of state parameter is 



„ I dark matter 
energy transfer < , , 

"■^ dark energy 



dark energy 
dark matter 



Pip 

The modified Klein-Gordon equation follows from Eq. ([4]): 



dip if 

When wc include baryons {pb) and radiation {pr)^ the remaining evolution equations are 



Pb 

Pr 

H 



-4.Hpr , 
2 



Pb 



r^Pr 



subject to the Friedman constraint, 



K^p 

3IP 



(5) 
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We can define effective equation of state parameters for the dark sector, which describe the equivalent uncoupled 
model in the background: pc + 3H{1 + Wc,cs)Pc = 0, Pip + 3H{1 + w^^cs)Pip = 0. By Eqs. ^ and dH), 



Wc,cS ■ 



iHp, 



3Hp^ 



It follows that 



> 
< 



Wc.cS > dark matter redshifts faster than a 

w<p,cff < ui^ dark energy has more accelerating power 

Wc,off < dark matter redshifts slower than a^'^ 

Wip^cS > Wip dark energy has less accelerating power 



(12) 

(13) 
(14) 



When Q > it is possible that w^p^cS < — 1 (see 01 for specific examples). This means that the coupled quintessence 
behaves like a "phantom" uncoupled model - but without any negative kinetic energies. 
Equations ([5|)- pil)) are an autonomous system of the form 



X = f (x) , 



(15) 



and the critical points satisfy f (x*) = 0. In order to study the stability of the critical points, we expand about them, 
X = X* + u, and Eq. (fTSj) yields 



Here g(x)/||x|| ^ as x ^ x*, and 



u-f'(x,)u + g(x) 



dx 



-(X*): 



(16) 



(17) 



is a constant non-singular matrix, whose eigenvalues encode the behaviour of the dynamical system near the critical 
point. 

If a component of f can be written as a fraction u(x)/t;(x), then a critical point requires the vanishing of the 
numerator, u(x,) — 0. If the denominator also vanishes at the critical point, w(x,) — 0, then care is needed in 
obtaining the eigenvalues of the linearized system (fT5|) . Strictly speaking, the fraction ■u(x,)/u(x^) may not be well 
defined. However, it is still possible to obtain analytical results via analysis of the behavior of the eigenvalues of f in 
the limit u(x) 0. 
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II. MODELS OF THE DARK SECTOR COUPLING 



There is as yet no basis in fundamental theory for a specific couphng in the dark sector, and therefore any couphng 
model will necessarily be phenomenological, although some models will have a more physical justification than others. 
Various models of energy exchange have been considered. Some of these are simple functional ansatzes, such as 
Q (X a". However these models are incomplete: they cannot be thoroughly tested against observations, since one has 
no idea what the perturbation of Q should be. 

A satisfactory model requires at least that Q should be expressed in terms of the energy densities and other covariant 
quantities. Two simple examples represent two of the main types of model: 

(I) Q = y/2/3nPp,ip, (18) 

(II) Q = aHp, , (19) 

where p and a are dimensionless constants whose sign determines the direction of energy transfer, according to Eq. ([5|): 

^ f > ; r f dark matter dark energy 

a,p < ^„ energy transfer < , , j i (20) 

I < 1 dark energy dark matter ^ ' 

For model (II). the case a > corresponds to the decay of dark matter into dark energy, with decay rate aH. 

Coupling (I) may be motivated within the context of scalar-tensor theories [1, 0, [6| . Generalizations of this model 
allow for P = P{(p) and more general forms of V{ip) (see, e.g., Refs. HB). Couplings which generalize or are closely 
related to model (II) have been considered as well, see, e.g., Refs. [3| for Q/H — UcPc + oi^Px and Ref. [l^l for 

For simplicity, we neglect the baryons (which are not coupled to dark energy), and we neglect radiation (since we 
are mainly interested in the late universe). The Friedman constraint (jlip becomes 

n^ + n^ = i, (21) 

and the total equation of state parameter is given by 

wtot := = — — — = w^ri^ and ptot + iH[l + Wtot)/Otot = . (22) 

Ptot Pip + Pc 

The condition for acceleration is Wtot < —1/3. 

We introduce the same dimensionless variables x,y as in the uncoupled case 0, where 

and Eq. (PT|) implies that 

< f],^ = + < 1 . (24) 
In the new variables, the equation of state parameters are 

2 

'^v = — ^ J ""^tot = - . (25) 

At a critical point (a;*,?/*), it follows that 

a(<)oct2/3(i+-^-yJ) and h > \i xl ~ yl < (26) 
The Hubble evolution equation may be written as 

§2=4^^ + -'-y')- (27) 

The energy balance equations ([3]) and ^ for coupling models (I) and (II) are independent of H when expressed in 
the variables x{N) and y{N), where = Ina. Thus Eq. (P7|) is not needed for these coupling models, and the phase 
space is two-dimensional (a;, y) space. 
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Coupling Model (I): Q = k. Pp^ip 



The autonomous system is 



(28) 
(29) 



where a prime denotes d/dN. The critical points and the conditions for stabihty, acceleration (wtot* < ~l/3) and 
physical existence, are summarized in Table HI where for convenience we have introduced the parameters 



6.= A-^a B _ V6M±(&^- 3/2)^/2 39/8)] 
3 ^ ' ^ 462 + 3/4 



(30) 



Table 
work 



J combines all the possible cases (A and (3 negative and positive), and recovers the particular results of previous 
H, [iH in the case of pressure- free matter {wc = 0). 



Point 




y* 


Stable? 






Acceleration? 


Existence? 


A 


1 





/3 > - and A > ^ 


1 


1 


no 


all A, P 


B 


-1 





3 

/3 < -- and A < -^6 


1 


1 


no 


all A, P 


C 


A 


M 


< A^ < 6 and 

A^ - 4=-^-3 < 0; 
A = 


1 




A^ < 2 


A^ < 6 


D 


26 


y^9-2\/6/36 


6^ > - and: 
- 2 

Bsgn(6) < Pb < or 


9-^/36 
362 


36 


6 6 


6^ > - and 
- 2 

f(3-!,")<^i,<2^ 


E 


2/3 
3 





stable when it exists 
and is hyperbolic 


4/3^ 
9 


4/3^ 
9 


no 





TABLE I: The properties of the critical points for the coupling model (I). Here 6 and B± are defined in Eq. (|30p . 



The complicated stability conditions for critical point D were confirmed numerically, and the results are shown in 
Fig. [H This point is an accelerated scaling solution that allows for Eq. ([T]) to be satisfied, i.e., 

< ilc* , < 1 and Wtot* < "'^ ■ (31) 

Acceleration requires (3/h < — %/6/6. This includes positive and negative (i (provided that A is accordingly restricted), 
but in the uncoupled case, /3 = 0, acceleration is not possible for the scaling solution. 
Point E is also a scaling solution, but it is always decelerating. 
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FIG. 1: The region of stability (blue, shaded) in the {b, (3) parameter space, for critical point D in coupling model (I). 

Coupling Model (II): Q = aHp^ 

The autonomous system of equations in this case is 

= _Sx + X^y' + Ul + x'-r/) + a^-^^-^-y^, (32) 
2 2 2x 

y' = -X^xy+^y{l + x'-y'). (33) 

We summarize the critical points and their properties in Table |lll Our results cover all signs of A and a and agree 
with previous results [ll| for A < 0, a < 0. Two additional critical points, E and F, occur for a > 0. 



Point 


a::. 


y* 


Stable? 




Wtot* 


Acceleration? 


Existence? 


A 


1 





a > 3 and A > \/6 


1 


1 


no 


all A, a 


B 


-1 





a > 3 and A < —VE 


1 


1 


no 


all A, a 


C 


A 

71 




Q > A^ - 3 and A^ < 6 


1 


3 


A^ < 2 


A^ < 6 


D 


Q + 3 


\/(q + 3)2 - 2qA2 


see Fig. [2] 


(a + 3)^ -aA^ 


a 


a < -1 


< 9 and 
2a<("+5<a + 3 


\/6A 


a/6A 


3A2 


3 


E 







A > and 


Q 

3 


Q 

3 


no 


< Q < 3 






C < a < 3 










F 







A < — \/6 and 
C < a < 3 


Q 

3 


Q 

3 


no 


< Q < 3 



TABLE II: The properties of the critical points for the coupling model (II). Here C A^ — 3 + ^ \^{\^ + 6). 

Furthermore, we have performed numerical integrations to probe the complicated stability conditions for the critical 
point D. The results are summarized in Fig. O Point D again allows for accelerated scaling solutions that satisfy 
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FIG. 2: The region of stability (blue, shaded) in the (A, a) parameter space, for critical point D in coupling model (II). 

Eq. ID), provided that a < —1. In the uncoupled case a = 0, and for positive a (decay of dark matter into dark 
energy), acceleration is not possible for D. 

Points E and F are also scaling solutions, but they are always decelerating. 



III. COUPLING MODEL (III): Q = Fp^ 

Coupling model (I) was introduced via scalar-tensor theory Q , and was shown to produce accelerated scaling 
solutions. Although model (I) has a clear physical motivation, it is ruled out by observations Q: the accelerated 
scaling attractor is not connected to a matter era where structure grows in the standard way. Indeed, generalizations 
of (I) with P = P{ip) face the same problem 

Coupling model (II) is not based on a physical model of dark sector interactions, but is designed for mathematical 
simplicity. This model and its simple generalization 0], Q = aH{pc + px), are specifically designed to produce an 
accelerated scaling attractor. Model (II) and its generalization evade the problems that model (I) has with a non- 
standard matter era These models are useful for phenomenology, but it is difficult to see how they can emerge 
from a physical description of dark sector interactions. The rate of transfer, aH, is determined by the expansion rate 
of the universe - rather than by purely local quantities associated with particle/ field interactions. 

In order to avoid this problem, we follow [1J| and replace the non-local transfer rate of model (II) by a local rate F, 



(III) Q = Fp, 



(34) 



where we assume that F is constant. This form of Q is used in other contexts. When F > 0, the same Q is used in: 

(1) a simple model to describe the decay of dark matter into radiation 

(2) a simple model for the decay of a curvaton field into radiation [l5l |. 

(3) a special case of a model in which superheavy dark matter partT^es decay to a quintessence scalar field JE]. 
When F > 0, the energy transfer in Eq. ([M)) corresponds to the decay of dark matter into dark energy. Models 

with decay of dark matter to dark energy allow for the possibility that there is no dark energy field in the very early 
universe, and that dark energy "condenses" as a result of the slow decay of dark matter. Coupling model (III) includes 
also the case F < which describes a transfer of energy from dark energy to dark matter: 




dark matter decays — > dark energy 

energy transfer from dark energy — > dark matter 



(35) 
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It turns out that the resulting evolution equations do not allow a two-dimensional representation of this model, since 
we cannot eliminate H from the energy balance equations ^ and (g]), using the variables x{N),y{N). Equation ([77|) 
must therefore be incorporated into the dynamical system. We do this via a new variable z, chosen so as to maintain 
compactness of the phase space: 

(36) 



H + Ho 



Thus < ^ < 1, and the compactified phase space now corresponds to a cylinder of unit height and radius. We also 
re-scale to a dimensionless coupling constant: 



7=7r- (37) 



Then we arrive at the autonomous system 



x' = -3x + \^y'+^-x{l + x'~y')~j^^^^^-^, (38) 

y' = ~X^xy+^y{l + x^-y^), (39) 

z' = ^z{l-z){l + x'-y^). (40) 

Since the system is invariant under y —y, the phase space may be reduced to a unit semi-cylinder. 

In order to determine the critical points and their stability, we need to deal with the singularities at a; = 0, z = 1 in 
Eq. (j38p . We rewrite the right-hand side of Eq. (j38p with a common denominator, in the form a:' = u(x, y, z)/[x(z — 1)]. 
Of the total of seven points that give y' — — z' and u ~ 0, three of them also give z — 1. The points A, B, C, D 
of Table UTTl are the critical points for the early universe, since z ^ corresponds to iJ ^ oo. These are the same 
critical points as occur in the uncoupled case. The new points E, F, G with z = 1 are late-universe versions of A, B, 
C {z^l => H ^0). 

Our results are summarized in Table IIIII The stability behaviour shown in Table IIIII is based on the eigenvalues 
of the linearized matrix, Eq. (fT7|) . which are shown in Table HVl In fact, only the signs of the real and imaginary 
parts of the eigenvalues are of importance. Therefore, we consider the eigenvalue ±oo simply as a positive/ negative 
eigenvalue. 

Apart from the extra dimension in its phase space, model (III) differs from models (I) and (II) in one key aspect: 

• In models (I) and (II) , the new behaviour introduced by coupling includes the possibility of an accelerated scaling 
solution (point D), characterized by Eq. ([T|). In model (III), no such accelerated scaling solution is possible. 
Instead, the accelerated attractor introduced by the coupling is point G, which has 

n„^0, fl^^^l. (41) 

This is the same qualitative behaviour as the standard ACDM model. 

Our analytical results are supported by numerical integrations, giving a consistent picture of the dynamical prop- 
erties of coupling model (III). Illustrative examples are shown in Figs. [3]and[31 Trajectories for A = 1 are shown in 
Fig- El showing the critical point G in Table [1111 In Fig. [H we plot trajectories when A = 4, and the stable node that 
corresponds to the point E is apparent. 

An interesting feature of the trajectories before reaching the global attractor (E, F or G, depending on the values 
of A and 7), is that they seem to be focused near a point that is vertically above the early- universe critical point D, 
and close to the z = 1 surface. Analysis of the system shows that the coordinates and derivatives of this (non-critical) 
point are 



x'^O, y'^Q, z' = 0(7)- (43) 

In the limit of zero coupling, we have z' ^ 0, and this point collapses to the z — Q critical point D. The deviation 
of this point from being critical is 0(7). When 7 is small, as in the plots, this explains the presence of focusing of 
trajectories in the numerics. 



Point 




y. 


2, 


Stable? 






Acceleration? 


Existence? 


A 


1 








saddle node for A > \/6 
unstable node for A < 


1 


1 


no 


all A, 7 


B 


-1 








unstable node for A > 
saddle node for A < — \/6 


1 


1 


no 


all A, 7 


C 


A 


V 







saddle node 


1 




A^ < 2 


A^ < 6 


D 


VE 

2A 


2A 





2 24 

saddle node for 3 < A < — 

2 24"^ 

saddle focus for A > — 
7 


3 

A2 





no 


A^ > 3 


E 


1 





1 


stable node for A > \/6 and 7 > 
saddle node for A > y/Q and 7 < 
saddle node for A < \/& and all 7 


1 


1 


no 


all A, 7 


F 


-1 





1 


saddle node for A > — \/6 and all 7 
stable node for A < —a/6 and 7 > 
saddle node for A < — \/6 and 7 < 


1 


1 


no 


all A, 7 


G 


A 

V6 


V 


''4 


1 


stable node for 7 > 
saddle node for 7 < 


1 




A^ <2 


A^ <6 



TABLE III: The properties of the critical points for the coupling model (III). 



Point 




y* 


z. 


Eigenvalues 


A 


1 








3; 3; 3-^^ 


B 


-1 








3; 3; 3+^^ 


C 


A 









D 


\/6 
2A 


2A 





1 -^(A±V24-7A2) 


E 


1 





1 


3; 3 2 ; sgn(7)oo 


F 


-1 





1 


-3; 3+— 2" ; -sgn(7)oo 


G 


A 

V6 


(■4)"' 


1 


-y; Y -3; -sgn(7)oo 



TABLE IV: Critical points and associated eigenvalues for coupling model (III). The 00 appears due to the limit z 
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IV. CONCLUSIONS 



We considered the background dynamics of a universe dominated by dark energy (in the form of exponential 
quintessence) and cold dark matter, where there is energy exchange in the dark sector, as in Eqs. ^ and 

Pc + 3Hpc = -Q = -[Pip + 3H{1 + w^)p^] . 
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For the previously introduced forms of Q, given in Eqs. ((T8)) and (fTQ]) . 

(I): Q= ^/2/3k(3p,^, (II): Q^aHp,, 

the phase space remains two-dimensional, as in the uncoupled case Q = 0. We found the properties of the critical 
points for all signs of A, a, (3. The results are summarized in Tables HI and [Til and slightly extend previous work P.[6l.[Tl| 
in the case of pressure-free matter (wc = 0). In both models, critical point D is the cosmologically relevant point, 
because for nonzero coupling it includes accelerated scaling attractor solutions, 

< ilc* , < 1 and wtot* < ~ 3 ■ 

The stability behaviour of critical point D was investigated numerically, and is shown in Fig. [1] for model (I) and 
Fig. m for model (II). 

Our main results arc for a new coupling model [l3|, defined in Eq. ([M)) . 

(Ill): Q = rpc. 

This has a similar form to model (II), but is more physical since the transfer rate F is determined only by local 
properties of the dark sector interaction at each event, and is not dependent on the universal expansion rate. When 
F > 0, this new model has the same form as simple models for the decay of dark matter particles to radiation [l4[, or 
to quintessence [l^ , and for the decay of the curvaton field into radiation ^15'| . 

Model (III) requires a three-dimensional phase space, since the Hubble rate cannot be eliminated from the equations 
for x',y'. This makes the dynamics more complicated than for models (I) and (II). In particular, a new set of late-time 
critical points arises in (III) , and considerable analytical effort is required to identify these points and determine their 
stability properties. Our results are summarized in Tables IIIII and IIVI We performed numerical integrations of the 
dynamical system in order to confirm the analytical results, and examples of these integrations are shown in Figs. [3] 
andm 

The cosmologically relevant critical point is G, which allows for an accelerated critical solution (when < 2). 
However, this is not a scaling solution, since 

f7„ = , n^, = i, 

which is similar to the asymptotic behaviour of the standard ACDM model. This accelerated critical solution is an 
attractor when 7 > 0, i.e., for the case when dark matter is decaying to dark energy. Note that in model (II), the 
decaying dark matter case, a > 0, does not lead to any accelerated attractor (see Table II). Model (HI) with F > 
produces an interesting class of models where dark matter decays to dark energy - so that the primordial universe 
may have no dark energy - and where this decay eventually leads to dark energy dominance, independent of initial 
conditions (since there is an attractor). Although such models do not solve the coincidence problem in the standard 
way, they may provide a new approach to the broader problem of explaining why dark energy dominates over dark 
matter only late in the universe's evolution. 

The background dynamics for coupling model (III) show new features not present in the previously investigated 
models (I) and (II). In order to confront this model with observations, the cosmological perturbations with a dark 
sector coupling of form (HI) need to be investigated (see Ref. p^). 
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